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A TRANSFORMATION METHOD FOR DERIVING, FROM A PHOTOGRAPH, 
POSITION AND HEADING OF A VEHICLE IN A PLANE 

Robert K. Sleeper and Eunice G. Smith 
Langley Research Center 

SUMMARY 

Equations have been derived that transform perspectively viewed planar surface 
coordinates, as seen in a photograph, into coordinates of the original plane surface. 

These transformation equations are developed in terms of nine geometric variables that 
define the photographic setup and are redefined in terms of eight parameters. The 
parameters are then treated as independent quantities that fully characterize the trans- 
formation and are expressed directly in terms of the four corner coordinates of a refer- 
ence rectangle in the object plane and their coordinates as seen in a photograph. 

Vehicle position is determined by transforming the perspectively viewed coordinate 
position of a representative vehicle target into runway coordinates. Vehicle heading is 
determined from the runway coordinates of two vehicle target points. When the targets 
are elevated above the plane of the reference grid, the computation of the heading angle 
is unaffected; however, the computation of the target position may require adjustment of 
two parameters. Methods are given for adjusting the parameters for elevation and an 
example is included for both nonelevated and elevated target conditions. 

INTRODUCTION 

Many situations exist which require monitoring the take-off, landing, and taxi behav- 
ior of aircraft or the ground handling performance of surface vehicles without special 
vehicle instrumentation and often under adverse operating conditions. A fundamental 
source of this type of information is the planar trajectory of the vehicle being examined. 

In the past, the trajectories have been straight lines and such devices as a trailing wheel 
or a phototheodolite have been adequate for providing the necessary trajectory data. 
Recently, however, interest has been generated in measuring vehicle behavior where the 
vehicle undergoes lateral displacements and/or changes in its heading as attributed, for 
example, to operations under cross-wind conditions. Since these trajectories do not lend 
themselves to the straight-line measurement techniques, a more general measurement 
technique is needed. Motion pictures in stereo were considered, but to the authors’ 
knowledge derivation of position data from still stereo motion-picture frame pairs has 



been unsuccessful because of synchronization problems of the twin views. Single frame 
motion pictures of a vehicle trajectory contain most of the positional data that the stereo 
process provides; however, no method is known to exist for reducing such data. The 
development of such a method was the objective of this report. 

This paper develops and applies a transformation technique to the trajectory deter- 
mination problem. The technique transforms coordinates in a single perspective view of 
a plane into coordinates of the original plane. The transformation uses parameters that 
are functions of several geometric variables which could, with difficulty, be measured. 
However, if a known rectangular reference grid is present in the test plane, the method 
of this report allows easy calculation of the parameters directly from the photograph. 

The requirement of a grid generally is not severe since such a grid is already present 
on many test surfaces as, for example, the expansion joints that exist on most concrete 
runways. For surfaces where no grid is present, one can be provided by the simple 
installation of markers, raised to the elevation of the test plane if necessary, at each 
corner of a measiured rectangle. 

Transformation equations are derived, a technique for determining values of the 
parameters is presented and, by using two reference target points on the vehicle, the 
position and heading are computed. For those cases where the targets are not in the 
plane of the grid, it may be necessary to adjust the parameters. Two example applica- 
tions of the method have been included. 


SYMBOLS 

Measurements were made in U.S. Customary Units and values have been converted 
to SI Units. 

Cj jth tilted perspective projection coefficient (see eqs. (12) to (17)) 

f distance from lens node to film plane 

H height of lens node above object plane 

AH elevation of target above runway 

i reference rectangle corner index, 1 to 4 

j tilted perspective projection coefficient index, 1 to 6 

k variable index, 1 to 4 


2 


I 


N 

O 


o 

P 

P 

Uk 


(X,Y) 

(Xi,Yi) 


(Xl.Yl) 

(Xn,Yn) 


(Xo>Yo) 


(Xr^Yr) 

(XsX) 


AX 

(x,y) 

(^L»yL) 


(Xn.yn) 


lens node (figs. 1 and 2) 

foot of perpendicular to object plane and passing through lens node 
(figs. 1 and 2) 

projection of point O on image plane (figs. 1 and 2) 

projection of principal point of image plane on object plane (figs. 1 and 2) 
principal point of image plane (figs. 1 and 2) 
kth composite variable (eqs. (B9) to (B12)) 
coordinates in object plane 

coordinates of ith corner of reference rectangle in plane of reference 
rectangle 

coordinates of left target in object plane 
coordinates of reference point in object plane 
coordinates of lens node in object plane 
coordinates of right target in object plane 

coordinates of origin of image plane of an original projection system in 
object plane of a secondary projection system (see eqs. (B3) and (B4)) 

X-component of elevation error 

coordinates in image plane 

coordinates of ith corner of distorted reference rectangle in image plane 
coordinates of left target in image plane 
coordinates of reference point in image plane 
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(xp,ypj coordinates of principal point in image plane 
(^R’^r) coordinates of right target in image plane 
AY Y-component of elevation error 

d rotation angle between image and object planes of perspective projections 

6 angle of tilt (figs. 1 and 2) 

$ coordinate rotation angle between Y -axis and principal plane (fig. 1 

and eq. (23)) 

(f) coordinate rotation angle between y-axis and principal plane (fig. 1 

and eq. (18)) 

\p heading angle 

A bar above a variable denotes principal coordinates in the image or object planes 
and a bar over a coefficient denotes an adjusted coefficient. Primes are used to denote 
quantities of different coordinate systems. 

TRANSFORMATION DERIVATIONS 

This report describes a method for determining the position and heading of an object 
on a flat surface from a single overhead photographic view of the surface. The method, 
when applied to motion pictures, offers a valuable tool for measuring data to describe the 
trajectory of an airplane on a runway. For such an application the camera may be tilted 
with respect to the runway and the vehicle is free to undergo lateral displacements and/or 
changes in heading. To meet the needs of the method, the only requirements are a rectan- 
gular grid marked on the runway and designated target points on the vehicle. The method 
applies a transformation for determining the positions of the target points relative to the 
runway grid from the motion -picture frames without iteration. The following sections 
present the mathematical development of the transformation where the target points are 
assumed to be in the plane of the reference grid, the formulas developed for calculating 
the parameters of the transformation from coordinates in the film and runway coordinate 
systems, and finally, an extension to the basic method to include the condition where the 
target points move in a parallel plane above the runway. 
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Development of Transformation Equations 

A fundamental premise of the technique is that coordinates of points of the object 
plane and the image plane of a photograph are related by a tilted perspective projection 
and coordinate rotations in the image and object planes as illustrated in figure 1. Both 
references 1 and 2, for example, provide a discussion of the perspective projection. In 
figure 1, N denotes a simple lens node and ON and pN are defined to be normal to 
the object and image planes, respectively. The point p on the image plane is called the 
principal point. Point P is the projection of the principal point through the lens node onto 
the object plane and point O is the foot of the perpendicular from the lens node to the 
object plane. The plane formed by the lines ON and pN is called the principal plane. 
Line pNP is called the optical axis; Q is the true angle of tilt of the lens optical axis 
with respect to the vertical; H denotes the lens node height above the runway; and f 
denotes a lens -to -film distance. Also shown are angles of coordinate rotations 
and 0 in the object and image planes, respectively. The rotations aline arbitrary coor- 
dinate (without bars) systems with the principal plane (with bars) coordinate systems. 

The coordinates with bars are called principal coordinates. 

The origins of the coordinate systems shown in the figure have been positioned at 
points O and p for clarity. However, the same axes designations are used in the sub- 
sequent derivation for arbitrarily translated origins. In the translated systems, the 
coordinates of point p are denoted Xp,yp and the coordinates of point O are denoted 
Xo,Yo. 

With the use of figures 1 and 2, the following relationships can be written between 
coordinates in the image and object planes: 


X = X cos 0 + y sin 0 (1) 

y = -X sin 0 + y cos 0 (2) 

X = X cos $ - Y sin (3) 

Y = X sin «> + Y cos (4) 

X - xp X - Xo 


f cos 0 - (y - yp) sin 0 H 

y - yp (y - Yq) cos 0 - H sin 0 
f (y - Yq) sin 0 + H cos 0 
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where 


Xp = Xp cos 0 + Yp sin <p 
Yp = -Xp sin <^ + Yp cos (p 
Xq = Xq cos - Yq sin # 


Yq = Xq sin $ + Yq cos ^ 

For the special case when the tilt angle 9 is zero, the image and object planes are 
parallel and the transformation reduces to a simple linear enlargement relationship. A 
method for deriving object plane data from image plane information for such a condition 
is straightforward; however, only the more general perspective transformations are 
treated herein. 

Equations (5) and (6) may be written as 


or 


X = Xo 


+ 


- ^p) 

(f cot 0 + yp) - y 





(7) 

( 8 ) 

(9) 


These transformation equations are seen to be characterized by nine independent quan- 
tities 0, 9, f, H, Xp, yp, Xq, and Yq which henceforth will be called per- 

spective variables. 

Let the transformation equations (7) and (9) be redefined in terms of constant coeffi- 
cients as 


X , C5 + 

Cl - y 

(10) 

^ Cl - y 

(11) 
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where 


Cj = f cot 0 + yp 

( 12 ) 

^ sin 6 

(13) 

C — ^ 

sin 6 

(14) 

0 

11 

(15) 

C5 = Xo 

(16) 

C 0 = Yq - H cot e 

(17) 


The six coefficients Cj, C 2 , Cg, C 4 , Cg, and C 0 are expressed here as 
simple functions of the perspective variables and determine the coordinate relationships 
of the perspective projection. In the next section, formulas are presented for computing 
these coefficients directly from coordinates measured in the image and object planes. 

Note that equations (12) to (17) define six coefficients in terms of seven perspective 
variables; thus, it is necessary to specify a perspective variable before all the variables 
that completely characterize the photographic setup can be ascertained from a picture. 

Solution for Transformation Parameters 

In this section the parameters characterizing the transformation equations (1) 
to (4), (10), and ( 11 ) are expressed explicitly in terms of a set of four point coordinates 
measured in both the image and object planes. The parameters consist of two angles 
associated with the rotational coordinate transformations and six coefficients which are 
related to the perspective projection. 

Solving for the transformation parameters requires that a reference rectangle in 
the object plane with sides parallel to an arbitrarily assigned X,Y axis system be 
visible in or be otherwise derivable from the image plane as illustrated in figure 3. Let 
the corners of the reference rectangle be assigned numbers as shown in the figure in a 
counterclockwise direction with the sides 1-4 and 2-3 parallel to the X-axis. The corner 
coordinates, so identified in both the image and object planes, are the only inputs required 
for the parameter solution. The details of how the angles and coefficients are determined 
from the set of input coordinate pairs are presented in appendix A; however, their derived 
expressions are given here with comments relative to the development. 
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Image plane coordinate rotation angle 4>: 

iy _ b _ y ^ 

(p = tan"l 


% % 

1-^ 


(18) 




where 


^a = 


_ (xg - xi)(xg - X4 ) (y 4 - Yi) + xi(x 3 - X 4 )(y 2 - Yi) - X 4 (x 2 - xi)(y 3 - Y4) 
(^3 - ''4)(y2 - yi) - (^2 - ^i)(y3 - y4) 


(19) 


ya = 


_ ~(y2 - yi)(y3 - y4)(^4 - ^i) - yi(y3 - y4)(^2 - - yi)(''3_- ^4) 


(^3 - ^ 4 )( y 2 - yi ) - (^2 - ^ i )( y 3 - y 4 ) 


( 20 ) 


Xb = 


_ (X 4 - X;^)(Xg - X2)(Y2 - yi) + Xi(xg - X2)(Y4 ' yi) - X2(x4 - Xi)(y 3 - 


(X3 - ^2)(y4 - yi) - (^4 - ^i)(y3 - y2) 


( 21 ) 


-(y4 - yi)(y3 - y2)(^2 - ^1) - y i(y3 - y2)(^4 : ^1) + y2(y4 -_^i),(''3.- ^2) 
(X 3 - x2)(y 4 - yi) - '(^4 - ^i)(y3 - y 2 ) 


(22) 


and the numerical subscripts refer to the corners of the distorted reference rectangle of 
figure 3(b). The denominator of equations (21) and (22) can vanish without inhibiting a 
solution for the angle 4>. This vanishing denominator, however, implies that the two 
sides of the quadrilateral in the image plane, 1-4 and 2-3, are parallel. Note, however, 
that a similar vanishing of the denominator of equations (19) and (20) cannot exist since 
it would implY that sides 1-2 and 3-4 are parallel. Once <p has been determined, equa- 
tions (1) and (2) maY be applied to anY image plane x,y coordinates to Yield principal 
x,y coordinates which will be used in the remainder of the derivation of the transforma- 
tion parameters. 


Object plane coordinate rotation angle 


= tan 


-1 


' Yl - Y2) 

X3 - ^21 


(y4 - yi)(y3 - ^2) 


| 1 / 2 | 


( y 4 - y 3 )(yi - y 2 ) 


- Y2' 

sgni IT 

W - X2y 


(23) 


Once $ has been found, equations (12) and (13) may be applied to any object plane 
X,Y coordinates to yield principal X,Y coordinates. Equations for determining the 


8 


perspective projection coefficients are presented, A general set of equations is given 
first; however, because the equations are indeterminant when ^ equals zero, a separate 
set is included for this case. Subsequent references will be made only to the general 
equations. 

Perspective projection coefficients 0): When the object plane rotation angle ^ is 
not zero, the perspective projection coefficients are solved sequentially in the order of 
the numerical subscripts as follows: 


^ _ y3(yi - y2)(X3 - X 2 ) tan ^ - Yi(y 3 - y 2 )(Yi - Yg) 
^ (yi - y 2 ) (X 3 - X 2 ) tan ^ - (y 3 - y 2 )(Yi - Yg) 


(24) 


IX 

1 

(ys - y2) 

(Ci - yi) - ( 

^2 - Xs)(y2 - yi)(Ci - y3) 

1 

(^1 - ^ 2 ) 

ys - (^1 - ^3 

j)y2 + (X2 - X3)yj 


K - Yl){ 

^1 - yi> 

(^1 ■ y2) 

C2(y2 - yi 

) 


(25) 


(26) 


C2[ci(x2 - xi) + (xiy2 - X2yi)] + (x ^ - X2)(Ci - y{){c^ - y2) 

^2(y2 - yi) 


(27) 


^ _ C2(x2 - xi) - Ci(x2 - Xj) + X2Y2 - XiYi 

*^5 — — 

y 2 - yi 


(28) 


^ Ci(Yi-Y2)-(Yiyi - Yjya) 
^6 — — 
y2 - yi 


(29) 


Perspective projection coefficients (4> = 0): When the object plane rotation angle is 
zero, the perspective projection coefficients are solved sequentially, not in the order of 
the numerical subscripts given below, but in the order C 4 , C 5 , C 2 , C 3 , € 5 , and Cj 
where 


Ci = 


yi(Yi - Ce) + C 2 C 3 


Yl-Cg 


(30) 
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( 31 ) 


(yi - y2)(Xi - C5) 


X2 -xi 


C3 


^ (Yi - Y2 )(xi - C4)(x2 - C4) 
(Xi - C5)(X2 - %) 


(32) 


^ ^ (x2 - X3)(x4Xi - X1X4) - (x i - X4)(xgX2 - X2X3) 

^ ' (^2 - ^ 3 )(Xl - X4) - (xr-'x4)(X2 - X3) 


(33) 


(Xi - X4)(x3X2 - X2X3) - (X2 - Xj) (X4X1 - X1X4) 

5 - (X2 - X3)(Xi - X4)- - (xj - X4)(X2 - X3) 


(34) 


Cr = 


Yi(xi - C4) - C3(Xi - C5) 


XI - C4 


(35) 


Adjustment of Parameters for Target Elevation 

In the development of the transformation, it has been assumed that the target is in 
the same plane as that containing the reference rectangle. However, for many applica- 
tions the target will move in a plane parallel to and above this "reference plane" and the 
previously derived transformation is no longer strictly valid. For such conditions, it 
may be convenient to think of the object plane raised above the reference plane. A good 
example of an elevated target is the case where the target is a point on the wing of an air- 
plane on a runway. One recourse is to eliminate the elevation by raising the reference 
corner markers to the height of the target. However, for those cases where it is impos- 
sible or impractical to obtain coplanar conditions, an adjustment of the transformation is 
required. 

Examine first the nature of a target position error when the target is elevated above 
the runway plane that contains the grid. Let (X + AX) and (Y + AY) be the apparent coor- 
dinates of an elevated target as determined by the transformation and X and Y be the 
coordinates of the same target orthogonally projected onto the runway. The difference 
between the projected and apparent positions is designated the error due to an elevated 
target. Thus, for a lens node position (Xq,Yq) and lens node and target heights above the 
runway denoted by H and AH, respectively, the error components are given by 
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(36) 


AX = M(x -Xo) 

AY = ^(Y-Yo) (37) 

and the 

Total elevation error = [(AX)^ + (AY)3j ^ - Xq^^ “ ^o)^J ^ 


which is independent of the orientation of the object plane coordinate system. As one 
would expect, these relations reveal that the elevation error varies with the ratio of the 
target elevation to the camera height and the distance of the target from the lens foot- 
print coordinates ^Xq,Yq). In the subsequent adjustment procedures, the effective eleva- 
tion is reduced by analytically raising the grid plane to that of the target. 

An interesting and important fact may also be noted with regard to vehicle heading. 
Let the heading angle be calculated by the expression 


xjy = tan-1 


Yr - Yl 

Xr -Xl 


(39) 


where the subscripts refer to left and right points. By using equations (3), (4), (10), 
and (11), equation (39) can be written as 


x{/ = tan“^ 


-|pl(^R ;/^l) + C 4 (yL - yn) + - X^yj sin + C^{y^ - y^) cos ^ 

|pl(^R - + ^4(yL - yR) + ^lFr - %yL] cos ^ + C 3 (yR - 5 Fl) sin 


(40) 


Suppose that the reference rectangle, assumed to be on the runway, was raised to the 
elevation of the target points and the transformation was recomputed. It is apparent 
from equations (12) to (17) that the only coefficients that would be affected by the change 
in height H would be C 2 and € 5 . Since neither of these coefficients appears in 
equation (40) , it follows that target elevation introduces no error in the heading angle 
calculation. 

The remainder of the section gives three rigorous methods of compensating for 
elevation error by adjusting the transformation parameters. Each of these methods 
requires that the transformation parameters for the runway be solved by using the trans- 
formation as given and, with additional information which includes the elevation, two of 
the coefficients adjusted. 
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Lens node height H known. - In principle, the adjustment can always be made if 
the elevation AH and one of the perspective variables are provided. The height H of 
the lens node above the reference plane is, perhaps, the least critical variable and the 
easiest to measure. Since the height required by the transformation is now (H - AH), 

C 2 and Cg must be adjusted so that equation (13) becomes 


- H - AH 

Lo = : 

Sin 6 


H 


sin 8 


1 - 


H 


- r II 

-C211 --g- 


and equation (17) becomes 

Cg = Yq - (H - AH) cot 8 


(41) 


= (Yq - H cot 0 ) + AH 


1 

1/2 

lc2j . 



h/C2 




(42) 


where the bar above the coefficient denotes an adjusted coefficient. 

Point in the elevated plane known.- When the coordinate pairs of a single reference 
point on the elevated target plane are known, for example, ^Xj^Yn) (^n>^n)> another 
method for adjusting the transformation can be employed. The angles 0 and and 
the coefficients to €5 are first determined from the transformation coefficient 

equations (18) and (23) to (29), based on the runway grid. The coordinates of the ele- 
vated reference point may then be inserted into equations (1) to (4) and the following 
manipulated forms of equations ( 10 ) and ( 11 ) may be used to compute adjusted values 
for C 2 and Cg: 

(Xn - C5)(Ci - yj 


C 2 = 


Xn - C 4 


and 


^ ^ C 3 (Xn - C 5 ) 

6 Xn - C 4 


(43) 


(44) 
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Principal point position known .- Cameras typically used in photogrammetry have 
the location of the principal point (xp,yp^ indicated by fiducial marks on the edge of the 
film. Knowledge of the location of this point permits a solution for all the perspective 
variables without additional information; when the elevation AH is supplied, the trans- 
formation can be adjusted for the elevation change in effective camera height. For this 
condition the transformation parameters are first computed from the transformation 
coefficient equations (18) and (23) to (29) based on the runway grid. Next, the coordi- 
nates of the principal point are converted into principal coordinates by using equations (1) 
and (2) to yield (jcp,yp). If sufficient photographic fidelity is present, the new value of 5Ep 
should agree with the value for C4. To adjust the coefficients C2 and €5 for eleva- 
tion, equations (12) and (14) may be combined to yield 


cos Q = 


^1 -yp 
C3 


(45) 


Equation (45) may be introduced into equations for the adjusted coefficients C2 and C0 
as shown: 


C2 


H - AH 
sin Q 


and 


= C2 


AH 
sin 0 


= C2 



AH 



Ce = Yq - (H - AH) cot 0 


= C0 + AH cot 0 

(Ci - yp) AH 


- C0 + 


fr2 ic V \2ll/2 

ys - yi- ypj _ 


(46) 


(47) 


Equations (46) and (47), with known values of yp and AH inserted, replace values of 
C2 and C0 derived from the runway grid. 


Effect of Image Plane Projection on Usage 

A single projection relationship has been assumed in the derivation of the preced- 
ing equations; however, most image planes are derived from more than one projection. 
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For instance, pictures which have been enlarged in photographic processing or film that 
has been projected in some other manner, such as for still or motion-picture projection, 
have been derived from two perspective projections. To show whether the perspective 
projection relationship between the image and object planes still exists for a projected 
image plane, transformation equations have been formulated in appendix B for two con- 
secutively applied tilted projections. The parameters of the transformation were found 
to change with the additional projection but the form of the equations was unchanged. It 
was concluded that the perspective projection relationship is retained for additional pro- 
jections and, hence, the transformation method is valid even when the photographic data 
have been altered by projection techniques. 

APPLICATIONS AND DISCUSSION 

Two applications of the transformation have been made to exercise and evaluate the 
method. One application had targets in the same plane as the reference grid to test the 
direct application of the transformation. The other application involved an elevated tar- 
get and required adjustment of the coefficients. 

Target and Grid Plane Coincident 

The main objective of this test was to examine the accuracy of the direct applica- 
tion of the transformation method to photographs obtained from different size cameras. 
The reference rectangle selected for this example was 2.743 m wide and 4,572 m long 
and targets were designated as shown in figure 4, one of the data photographs. The 
camera lens node position for each photograph was approximately 4.3 m above the sur- 
face with approximate coordinates of (-2.7 m,-3.2 m) measured from the lower corner 
of the reference grid. The camera tilt was approximately 54°. Four general inventory 
still cameras of types standard for their format size were employed to represent four 
different camera format sizes and lens-to-film distances. Photographs were made with 
each camera and, as required, enlarged to a common size of 20 cm by 25 cm. In the 
cases where enlargement was required, the image plane was then derived from a double 
projection as discussed in a previous section. 

The transformation equations were programed on a desktop calculator equipped 
with a digitizer to convert points on the photograph electronically into x,y coordinates 
to the nearest 0.25 mm. 

Results are presented in table I which lists the actual coordinates of the two tar- 
gets in the plane of the reference grid and the corresponding coordinates determined 
from the photographs by using the transformation directly. The maximum error is 
found to be 0.8 cm where the width of the grid markings was 1,3 cm and the errors are 
seen to vary little for different cameras. 
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Targets Elevated Above the Grid Plane 

The objective of this test was to measure from 16 -mm motion-picture film the 
position and heading of an aircraft landing-gear model on a runway when the vehicle 
targets are elevated above the runway. By use of a general inventory camera, motion 
pictures were taken of the model placed at nine different positions on the runway and 
oriented with heading attitudes of -30°, -15°, 0°, 15°, and 30°. Frames from this 
film are presented in figure 5 which shows the model with a heading attitude of -30° at 
each of the tested positions. Two circular targets 6.4 cm in diameter and elevated 
0.271 m above the runway are identified on the model. (The targets were spaced 1.219 m 
apart and attached to outriggers to facilitate heading measurements.) Grids in the cen- 
tral section of the runway had a 1.219-m lateral and longitudinal spacing and corner 
coordinates as specified by the grid in the figure. In this coordinate system, the motion- 
picture camera which was 3.246 m above the runway had an approximate location of 
(0.7 m,3.3 m). 

The steps in determining the model position and heading were straightforward. 

The motion pictures of the tests were first projected onto the digitizer platen; this pro- 
cess enlarged the frame size to approximately 0.37 m by 0.52 m and then the film move- 
ment was stopped. (Like the first example, this application necessarily incorporated a 
double projection.) The coordinates of the four grid corners (indicated with circles for 
each position in fig. 5) and two targets shown in the projected picture were read by the 
digitizer and input directly into the programmable desktop calculator. Transformation 
parameters were computed from both of the image and object plane grid corner coor- 
dinates by using equations (24) to (29); then, by introducing the known camera height and 
target elevation into equations (41) and (42), the parameters were adjusted for elevation. 
Target point coordinates were computed by inserting the digitizer target readings into 
the transformation equations. Finally, the heading angle was determined from the tar- 
get coordinates by using equation (39), and the model position was found by averagii^ 
the two coordinates. 

Test results were computed from digitizer readings made by an operator unfamiliar 
with the test and no attempt was made to enhance original readings. These results are 
presented as determined vehicle positions and headii^s in table II. When compared with 
their actual values, most position differences are less than 1 or 2 cm and most head- 
ing discrepancies are less than a degree. (The width of the grid lines was about 3 cm.) 
Inaccuracies found in the table could generally be attributed to the poor photographic 
clarity of the films; nevertheless, the results obtained were judged to be sufficiently 
accurate for evaluations of landing-gear performance for which the method was intended. 
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General Comments 

This photographic measurement technique has a number of advantages that should 
be emphasized. The technique utilizes basic equipment that is readily available, it does 
not interact with the vehicle being measured, and, because the technique does not require 
the vehicle to carry onboard equipment, it can be used equally well for full-scale vehicles 
and for models where weight and volume are limited. Furthermore, little advanced prep- 
aration of the model or the surface is required. 

Special emphasis should be placed on the use of the technique when no elevation 
adjustment is necessary. Not only is less computation required but, for such conditions, 
only the reference grid coordinates need to be designated, and all the measurements for 
computing the transformation parameters can readily be made from the film. Since all 
measured transformation data are now inherent in the film, there is more photographic 
freedom. For instance, photographs need not be taken from a camera whose position is 
fixed, but moving reference bases, such as those provided by aircraft, where speed and 
attitude may be varying or unknown may be used. Furthermore, photographs may be 
taken with cameras employing such devices as zoom lenses which vary the lens-to-film 
distance. When there is no need for elevation adjustment, the only requirement of the 
technique is that a known reference rectangular grid be in the object plane and the field 
of view of, or otherwise derivable from, the photograph. This requirement is not con- 
sidered to be severe since such a grid is already present on many surfaces as, for 
example, the expansion joints that exist on most concrete runways. For surfaces where 
no grid exists, one can be simply introduced with the addition of markers, elevated to the 
target height if necessary, at each corner of a measured rectangle. With special mar- 
kers this method may even be applied to measure the position of a vehicle on an unusual 
surface such as water. 

It should be noted, however, that errors may exist in the image plane due to photo- 
graphic distortion, digitizer resolution, and reading inaccuracy; and such errors are pro- 
portionally amplified by the height of the camera above the grid plane. Another type of 
error multiplication also exists which is not bounded. This error occurs for points near 
the horizon trace where a position difference in the object plane can be imperceptible in 
the image plane and is caused by the singularity of the transformation in that region. For 
this reason data measurement near the horizon of a photograph should be avoided. 

CONCLUDING REMARKS 

Equations have been derived that transform coordinates from a perspectively viewed 
planar surface into coordinates of the original plane surface. These transformation equa- 
tions are developed in terms of nine geometric variables that define the photographic setup 
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and are redefined in terms of eight parameters. The parameters are then treated as 
independent quantities that fully characterize the transformation and are expressed 
directly in terms of the foim corner coordinates of a reference rectangle in the object 
plane and their coordinates as seen in a photograph. 

Vehicle position is determined by transforming the perspectively viewed coordinate 
position of a representative vehicle target into runway coordinates. When the target is 
in the plane of the reference grid, only the coordinates of the grid and the four corner 
coordinates from the perspective view are needed to derive the transformation. However, 
if the planes of the target and grid are not coincident, it may be necessary to adjust two 
of the parameters of the transformation for the elevation difference. Three methods for 
adjusting the parameters for an elevated target are given. Vehicle heading, which is 
found by deter minir^ the runway coordinates of two selected points on the vehicle that 
have the same elevation, is shown to be independent of elevation difference. 

Two tilted arbitrary projections applied in sequence were found to be equivalent to 
a single projection, from which it was deduced that the transformation method is valid 
even when there are supplemental projections as are common in the photographic process 
and occur anytime that still or motion-picture film is projected. 

Several advantages of the technique are evident. Some of these advantages are 

(1) It measures all planar trajectory data without interacting with the vehicle being 
measured and carries no onboard equipment. 

(2) It requires little preparation of vehicle or surface for its application. 

(3) It is equally useful on model or full-scale tests. 

(4) It utilizes simple, inexpensive photographic equipment that is readily available. 

(5) It permits the use of photographic projection techniques for improved data 
analysis. 

(6) It requires only minimal computation capability. 

The example applications of the technique showed little variation in the results for 
cameras with different size formats, and the technique was judged to be potentially use- 
ful for ground vehicle performance analyses. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
March 11, 1976 
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SOLVING FOR THE TRANSFORMATION PARAMETERS 

In this appendix transformation parameters (p, C^, C2, C3, C4, C5, 

and C0 are derived in terms of a set of four corresponding image- and object-plane 
coordinate pairs for coordinates in the object plane positioned in a rectangular pattern. 
Figure 3 (a) shows a reference rectangle arbitrarily oriented in the object plane. An 
X,Y axis system is established with axes parallel to the sides of the rectangle as shown 
and corresponding corners of the reference rectangle are numerically designated in a 
counterclockwise order with sides 1-4 and 2-3 parallel to the X-axis of the object plane. 
Also shown is a rotated X,Y axis system with its ordinate parallel to the principal plane 
intersection which will be determined. The amount of clockwise rotation of the X,Y axis 
with respect to the object plane X,Y coordinate system is denoted by 

An image plane view of the reference rectangle of figure 3 (a) is shown in figure 3 (b). 
The coordinates of the four corner points of this quadrilateral may be referred to any con- 
venient x,y axis system. In the image plane it can be seen that, in general, sides of the 
quadrilateral exhibit vanishing points such as those at a and b. These points are on the 
horizons of a photograph and define a line called the horizon trace. 

Also shown in figvtre 3 (b) is an x,y axis system rotated counterclockwise by an 
angle <P with respect to the x,y axis system. This axis system has its x-axis par- 
allel to the horizon trace as shown; its y-axis must be parallel to the intersection of the 
image and principal planes; therefore, the angle <P specifies the orientation of this 
intersection in the image plane. 

The rotational transformations associated with the angles (p and 4 > convert the 
image and object plane coordinates to their respective principal coordinates. After the 
rotational transformations in the image and object planes have been determined, the coef- 
ficients of the perspective projection may be evaluated in terms of the known image and 
object plane coordinates converted to the principal coordinates. 

The sets of four coordinate pairs then are the only requirements for a solution to 
the transformation parameters. These sets consist of the corner coordinates ^XjjYjj, 
(X2,Y2), (X3,Y3), and (X4,Y4) of the reference rectai^le in the object plane, and the 
corner coordinates (xi,yi), (x2,y2)> (x3>y3)> and (x4,y4) of the distorted rectangle as seen 
perspectively in the image plane. The first parameters to be determined are the angles 
of the rotational transformations. 
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Determination of Image Plane Rotation Angle 

The image plane rotation angle <p is determined from the image plane coordinates 
by finding the vanishing points where sides of the quadrilateral intersect, by constructing 
the horizon trace line that passes through these points and by determining its slope. This 
slope corresponds directly to the angle <p. 

Coordinates of the vanishing points are found by solving for the points of intersec- 
tion of the lines that pass through adjacent corner points and are given by 


(X2 - Xi)(x3 - X4)(y4 - y^) + xi(x3 - X4)(y2 - Yi) - X4( xa - Xi)(y3 - y4) 

(^3 - '^4)(y2 - yi) - (^2 - ^i)(y3 - y4) 


(Al) 


-(y2 - yi)(y3 - - ^i) - - y4)(^2 - ^i)-*-y 4(y2 -yi)(^3 -^4) 

(^3 - ^4)(y2 - yi) - (^2 - ^i)(y3 - y4) 


(A2) 


N i^)(^3 - ^2)(y2 - yi) + ^1(^3 - '^2)(y4 - yi) - ^2(^4 - ^i)(y3 - y2) 

(X3 - X2)(y4 - yj) - (X 4 - Xj)(y3 - y 2 ) 


(A3) 


-(y4 - yi)(y3 - y2)(x2 - Xj) - yj(y3 - y2)(x4 - Xi) + y2 (y4 - yi)(x3 - X 2 ) 

(^3 - ^2)(y4 - yi) - (^4 - ^i)(y3 - y2) 


(A4) 


where (Xg^jy^) denotes the vanishing point associated with the intersection of the extended 
lines 1-2 and 3-4 and (xjjjyjj) is associated with lines 1-4 and 2-3 of figure 3(b). 

The equations for each vanishing point are shown to have a common denominator 
which is related to the difference between the slopes of opposite sides of the quadrilateral 
in the image plane. Thus, the denominator is zero and is said to have a vanishing point at 
infinity when opposite sides are parallel. However, for the specified orientation of the 
quadilateral only the sides 1-4 and 2-3 can be parallel (zero tilt excluded). Thus, only 
the denominators of equations (A3) and (A4) can vanish. 

The angle (p may be written as 


4> = tan“l 


yb-ya 

^b - ^a 


19 


I 



APPENDIX A 


However, an indeterminant form exists when the denominator of equations (A3) and (A4) 
vanishes and the alternate form 


0 = tan'l 


Xb ■ Xb \ 



(A5) 


which does not exhibit this characteristic will be more useful. 

Once <(> is determined, equations (1) and (2) may be applied to any image plane 
x,y coordinates to yield the principal x,y coordinates which are used in the derivation 
of the transformation parameters. 


Determination of Object Plane Rotation Angle 

For the solution of the object plane rotation angle <i>, it can be noted that the object 
plane coordinate system has been alined with the reference rectangle so that 


X4 - Xj = Xg - Xg (A6) 

Y4 - Yi = Yg - Y2 = 0 (A7) 

and 

Xj - Xg = X4 - Xg = 0 (A8) 

Yi - Y2 = Y4 - Yg (A9) 


Introducing a portion of relations (A7) and (A8) in equation (4) yields 


and 


or 


sin ^ 


Y3 - Y2 

X3-X2 


cos $ 


Y1-Y2 

Yl -Y2 


(X3 - Xj) (Y, - Yj) 


(AlO) 


(All) 


(A12) 
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where the quantity ^¥3 - Y2)y/^Yj - ¥2^ is an unknown ratio whose magnitude and sign will 
be determined. 

When equations (A6) to (A9) are introduced into equation (4), it is also found that 

¥4 - Yi = ¥3 - ¥2 (A13) 

and 

Yi-Y 2 =Y 4 -Y 3 (A14) 

From the identity (A13) and equation (11), the following relation may also be found: 


C2C3 C2C3 C2C3 C2C3 

Cl - y 4 Cj - yi Cl - y 3 Ci - y2 

or since the product C2C3 can never be zero, 

y4 - yi _ y3 - ^2 

(Cl - y4)(C, - 7,) ■ (Cl - y3)(Ci - yj) 

By using the identity (A 14 ) and equation ( 11 ) in a similar manner, 

yi - y2 . y4 - y3 

(Cl - yi)(Ci - y2) (Cl - y4)(Ci - yg) 


(A15) 


(A 16) 


(A17) 


Equations (A16) and (A17) may be combined to yield 


or 


(^1 - yif 

(y4 - yi)(yi - y2) 

\Cl - Y3/ 

(y 3 - y2)(y4 - y3) 


- yi 

Cl - yg 


(y4 - yi)(yi - y2), 


( y3 - y2 )( y4 - y3 ) 


1/2 


(A 18) 


(A19) 
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With the aid of equation (11) the unknown ratio of equation (A12) may now be expressed 
in terms of the ratio of equation (A19). Thus, 


Yi - Y2 = C2C3 


yi -y 2 

(Cl - yi)(Ci - yj) 


(A20) 


Y3 - Y2 = C2C3 


ys -y 2 


(Cj - y3)(Ci - yg) 


(A21) 


or 


^3 - Y 2 . (^1 - yi)(y3 - y2) 

Yi - Y2 (Ci - y3)(yi - y^) 


or when equation (A19) is inserted, 


Y 3 - Y2 

Y1-Y2 


= ± 


(y 4 - yi)(y3 - y2) 


(y4 - y3)(yi - y2) 


1/2 


(A22) 


(A23) 


Thus, upon substituting equation (A23) into equation (A12), the magnitude relationship 
of # may simply be expressed as 


tan-1. 


X3 -X2 


(y4 - yi)(y3 - y2) 


(y4 - y3)(yi - y2) 


1/2 


(A24) 


The sign of may be determined by examining the following sketches of the 
reference grid in the object plane and the grid in the corresponding principal image 
plane for different values of 4>. (See sketches (a) and (b).) These sketches show that 
parallel sides of the quadrilateral in the image plane correspond to zero object plane 
rotation (4> = 0) and the edges of the reference grid parallel to the X-axis of the object 
plane have a slope direction in the principal image plane commensurate with the object 
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Object plane 
Sketch (a) 



Image plane 
Sketch (b) 


plane rotation angle The sketches also indicate that for the specified corner desig- 
nation, the magnitude of the angle is less than 90°. Thus, 


$ = tan"l 


/Yl - Y2\ 

(y4 - yi)(y3 - y2) 

1/2 

CO 

1 

(y4 - y3)(yi - %) 



sgn 


yp - y2 
53 - ^2 


(A25) 


Once ^ has been specified, equations (3) and (4) may be used to convert object 
plane coordinates (X,Y) to the principal coordinates (x,y) for the solution of the coeffi- 
cients of the perspective projection which follows. A general solution for the perspec- 
tive projection coefficients has been derived first; however, it exhibits an indeterminant 
form when ^ is zero. A separate set of equations is derived for this condition. 
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Solution for Perspective Projection Coefficients Nonzero) 

Determination of Cj.- The coefficient Ci is the composite quantity ff cot 0 + yp) 
and, in the principal coordinate system of the image plane, it defines the horizon trace or 
limitii^ positive ordinate position. This coefficient may be found by combining equa- 
tions (A12) and (A22) to yield 


tan 


(^1 - Yz) 

(y3 - y2) 

(c 

(X3 - Xj) 

(yi - ^ 2 ) 



-Jl) 

-h) 


Thus, 


or 


y3(yi - y2)(^3 - ^2) tan 4> - yi (y3 - y2 )(Yi - Yg) 

(yi - y2)(X3 - X 2 ) tan «• - (yg - jF 2 )(Yi - Yg) 


Cl = 


y3(yi - y2)(^3 - ^ 2 ) - yi(y3 - y2)(^ 1 - ^ 2 ) cot ^ 
{fi - f2){^3 - ^ 2 ) - (y3 - y2)(^i - '^ 2 ) ^ 


(A26) 


(A 2 7) 


(A28) 


Determination of C2-- The coefficient C2 is the composite height function 
(H/sin 9) which is determined from three image-to-object plane coordinate pairs with 
the use of equation (10). Thus, for two corner coordinate pairs, for example, those 
identified by 1 and 2, 


Xi - Xg = Cg 


^1 - C4 
fl - fl 


X2j^\ 

Cl - y2/ 


and for points 2 and 3, 


X2 - Xg = 


- C4 

^\Ci - y2 


^3 - C4 ’ 

Cl - Y3J 


(A29) 


(A30) 


The coefficient C4 can be eliminated from equations (A29) and (A30) and, since Cj has 
been previously determined. 


C2 


{^1 ■ ^2) (^2 ■ yi)(^2 " ^3)(*^i ■ ys) ■ (^3 ■ y2)(^i ■ ^2)(ci - yi) 


(y3 - y2)(^iy2 - ^2yi) - (y2 - y^te - ^3y2) - Cil(' 


r/^ ^ 


^1 - ^2)(y3 - y2) - (^2 - ^3, 


!)(y2 - yi) 


(A31) 
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or after some further manipulation 


(Xi - X2)(y3 - Yz) 

(Cl - yj) - (: 

^2 - X3)(y 2 - yi)(ci - ¥ 3 ) 

(""I - '^2)y3 - ("^1 - ""a 

)y2 + (^2 - ^3)yi 


(A3 2) 


Determination of C3 .- The coefficient C3 is the composite lens -to -film distance 
function (f/sin 9) which can be determined by the application of equation (11) to two coor- 
dinate pairs, for example, those identified by 1 and 2. Thus, 


Yl 


- Y2 = C2C3 


yi - y2 

(Ci - yi)(Ci - Y2) 


or since Cj and C2 have been previously determined. 


(A33) 


(^2 - Yl)( 

, 0 

1 

1 

(^1 - y2) 

C2(y 2 - yi 

) 


(A34) 


Determination of C4.- The coefficient C4 is the coordinate Xp which is the 

abscissa position of the principal point in the principal coordinate system of the image 
plane. Since and C2 have been previously determined, equation (A29) may be 

readily rearranged to show that 


C 2[^l(^2 - xi) + (xjy2 - X2yi)] + (Xi - X2)(Ci - yi)(Ci - y2) 

C2(y2 - yi) 


(A3 5) 


Determination of C5.- The coefficient C5 is the coordinate Xq which is the 
abscissa position of the camera lens node in the principal coordinate system of the 
object plane. To find its value, equation (10) has been expressed in terms of the corner 
coordinate pairs identified by 1 and 2, and C4 has been eliminated from the relations 
to yield the expression 


C 


C2(^2 - ^1) - Cl(x2 - Xi) + X2y2 - XiYi 
Y2 - n 


(A36) 
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Determination of C0.- The coefficient Cg is a composite quantity (Yq - H cot 9^ 

which is the limiting negative ordinate position in the principal coordinate system of the 
object plane. It corresponds to a point at negative infinity for the principal coordinates 
in the image plane. In its solution, equation (11) may be expressed in terms of the coor- 
dinates identified by 1 and 2 and the coefficient C3 eliminated from the two equations 
to yield 


ce = 


Ci(Yi - Y2_) - (yi - Y2y2) 

Y2 - Yl 


(A3 7) 


Solution of Perspective Projection Coefficients (# Zero) 

When is zero, object plane coordinates (X,Y) are also principal object plane 
coordinates. Furthermore, certain principal image ordinate values are equal; that is, 

y2 = Y3 (A3 8) 

and 

y j = y^ (A39) 

By using these relations in equations (10) and (11), the perspective projection coefficients 
may be solved in the order that follows. 

Determination of C4 and C5.- Let equations (10) and (11) be expressed for the 
case when equals zero as 


Cl - y = 


C2C3 

Y - Cg 


Cl 


C2(x - C4) 
X- Cg 


(A40) 


(A41) 


Then, since yj = y4, equation (A41) yields 

Xi - C4 X4 - C4 

Xi - Cg " X4 - Cg 
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and similarly for Y 2 = Ys 

^2 - ^4 _ ^3 - ^4 
X2 - Cg Xg - Cg 

or 

C4(Xi - X4) - Cg(xi - X4) = X4X1 - X1X4 

and 

^4(^2 - X3) - Cs(x2 - X3) = X3X2 - X2^3 


These equations may be solved simultaneously to yield 


and 


C4 = 



- ^3) (^4X1 - X1X4) - (xi - X4)(xgX2 - XgXg) 

(^2 ■ ^ 3)(^1 ■ ^ 4 ) ■ (^1 " ^ 4)(^2 " ^ 3 ) 


-(Xi - X4)(xgX2 - X2X3; 

> + 

(X2 - Xg)(x4Xi - X1X4) 

(5^2 - - X4) 

-( 

1 

XI 

^2 - X3) 


(A42) 


(A43) 


Determination of € 2 *“ Let equation (A41) be applied to two corner coordinate 
pairs having different abscissa values in the object plane. Then, 

_ _ ^ f ^2 - C4 - C4\ 

- ^2 = = 2 ( 5 ^ - 

and since Xj = Xg and Cg is known 


C2 = 



(A44) 
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Determination of C3.- Eqxiations (A40) and (A41) may be used together on two 
corner coordinate pairs with different ordinates in the object plane to show that 


or 


Y1-Y2 



X2 - Cs\ 

X2 - C4/ 


- (Yi - Y2)(xi - C4)(x2 - C4) 

^ (^1 - C5)(x2 - h ) 


(A45) 


Determination of Cq.- Since C3 has been previously determined, equations (A40) 
and (A41) may be used together to solve for €5. Thus, 




Yi(xi - C4) - C3(Xi - C5) 


XI - C4 


(A46) 


Determination of Cj.- Since C2, C3, and C0 are now known, equation (A40) can 
be used to yield 


Cl 


yi(Yi - Ce) + C2C3 
Yl - Cg 


(A47) 


The eight transformation parameters that characterize the transformation have 
been solved in terms of four image and object plane coordinate pairs. A simple check of 
the transformation may now be performed by inserting the four image plane coordinates 
into the transformation equations and comparing the computed object plane coordinates 
with those used for the determination of the parameters. 
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EQUIVALENCE OF TWO SEQUENTIAL PERSPECTIVE 
PROJECTIONS TO A SINGLE PROJECTION 

In this appendix it is shown that two general perspective projections performed 
sequentially are equivalent to a single perspective projection. This relationship is 
proved by derivii^ the image-to-object plane relationship and showing that its equa- 
tions have the same form as those for a single perspective projection. 

Let equations (1) to (4), (10), and (11) be combined to give two explicit relations 
for the object plane coordinates in terms of the image plane coordinates. Thus, 


X = 


•|(CiC5 - C2C4] cos <i> + (CiC0 + C2C3) sin # 


Cj + (x sin (p - y cos 4 >) 

+ x|j^C2 cos (p + sin 0) cos # + C0 sin (p sin ^ 

+ y|^C2 sin (p - cos (p^ cos # - C0 cos <p sin 


(Bl) 


and 


Y = 


•^CjC0 + C2C3) cos 4* - ^CjCg - C2C^^ sin 


Cj + (x sin (p - y cos <p) 

+ xjc0 sin (p cos # - ^C2 cos 0 + C5 sin <p) sin ^ 
- y|c0 cos <p cos # + ^C2 sin (p - cos <p^ sin 


(B2) 


Let the image coordinates (x,y) be expressed in terms of new coordinates (X',Y') 
which are rotated by an angle a and translated by distances Xg and Y*g in a man- 
ner such that 


X = (x' - Xg) cos 0! + (y' - Yg) sin a 
y = -(x' - Xg) sin a + (y' - Y*g) cos a 


(B3) 

(B4) 
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Finally, let the coordinates (X*,Y*) be treated as the object plane of a second per- 
spective projection relationship distinguished by primes and given by 


X' = 


^.)({c'lC5 - c’aC^ cos + (c'lC^ + C’gC^) sin 


C|[ + (x' sin 4>' - y* cos 4>' 

+ x'jjc '2 cos <!>' + C 5 sin 4>') cos sin 4>' sin 

+ y* |^C 2 Sin <p' - Cg cos cos - Cg cos <t> sin 


(B5) 


and 


Y' 


^)[(ciC6 + C2C3) cos - (C1C5 - C’2Cy sin 


Cj + (x* sin <f) - y' cos (j>' 

+ x'|c0 sin 4 >' cos - (c'2 cos (p' + Cg sin </>') sin 

- y'|cg cos 0' cos + ^C2 sin 0' - Cg cos 0'j sin 


(B6) 


Expressions relating the original object plane coordinates (X,Y) and the image 
plane coordinates of the second perspective projection (x',y’) may be found by combin- 
ing equations (Bl) to (B6) and after considerable algebraic manipulation 


X 


= ([cslc’lUi - C^U2 - C^U4) + C2(c^U2 - ciUg - C;U4^ cos ^ 

+ [cg(CiUj - C4U2 - C3U4) + C'iC2C^ sin ^ + x' |jcg(Ui sin 0' + U2 cos 0*j 
- C2(U3 sin 0’ - U4 cos 0'^ cos $ +j(CgUj + C2C^ sin 0* + CgU2 cos 0^ sin ^ 
+ y'[[c5(u2 sin 0’ - Uj COS 0') + C 2 (U 3 cos <t> + U4 sin 0 '^ cos # 

+ |cgU2 sin 0' - (CgUj + C2C3) cos 0^ sin - C4U2 - C3U4) 

+ x'^Uj sin 0* + U2 cos 0') - y'(Uj cos 0* - U2 sin 0*)J 


(B7) 
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and 


Y = ^[c6(ciUi - C^U2 - C^U4) + CiCgCg] cos $ - [c5(ciUi - C^Ug - C^U4) 

+ ^2(^3^2 " ^1^3 " ^4^4^ ^ ■*■ ^2^3) ■'■ ^6^2 ^ 

- jc5(Uj sin cf>' + U2 cos <^') - C2(Ug sin 0* - U4 cos 0’^ sin ^ + y' ^ C0U2 sin 

" (^6^1 ■'■ ^2^3) ^ " [ps(^2 ~ ^1 ^') ■*■ ^2(^3 

+ U4 sin 0 ') sin - C4U2 - C3U4j + x'^Uj sin 0 ' + U2 cos 0 ') 

- y'(Ui cos 0 ' - U2 sin 0'^ ^ (Bi 


where 


Ui = Cl - Xg sin (0 + O') + Yg cos (0 + a) + C5 sin (4>' + 0 + a) 
- C0 cos (4>' + 0 + O') 

U2 = C'2 sin (#' + 0 + a) 


(B9) 


(BIO) 


113 = 04 + Xg cos (0 + a) + Yg sin (<p + a) - C*5 cos (4>' + 0 + a) 

- C0 sin (4>' + 0 + a) (Bll) 


U4 = C2 cos (#' + 0 + O') 


(B12) 


If, and only if, two consecutive perspective projections are to be equivalent to a single 
perspective projective, equations (B7) and (B8) must be of the form 


X = 


-((cicg - C^C^) cos *'• t (c’icf + C^Cj') sin «•' 


C'l + (x* sin 0" - y' cos 0' 

+ x'|(c2 cos 0" + C5 sin 0") cos 4>" + C0 sin 0” sin #'ij 
+ y'j^C'2 sin 0" - C5 cos 0") cos - €5 cos 0" sin 


(B13) 
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and 


Y = 


^|(CiC6 + CgCg) cos # - (C1C5 - C2C4) sin ^ 


Cj + (x’ sin <p" - y' cos 0 ") 

+ x' jcg sin 0" cos - (C2 cos 0" + Cg sin 0") sin 

- y'|c0 cos 0" cos + ^C'2 sin 0” + C5 cos 0") sin 


(B 14 ) 


where the double primed quantities denote equivalent transformation parameters of a 
sir^le perspective projection. 


It will now be shown that equations (B 7 ) and (B8) have the same form as equa- 
tions (B 13 ) and (B 14 ). To show this equivalence, it is convenient to first satisfy the 

denominator expressions by dividing the numerator and denominator expressions of 

\l/2 

Thus when like coefficients in the denomina- 


equations (B 7 ) and (B8) by (uf + u|) 
tor are compared, 


„ Ui sin 0 * + Vo cos 0' 
sin 0 = ■ 




2 \ 1/2 


and 


= 


Uj cos 0* - U2 sin 0* 


(uf 




or 


0” = tan-1 


Uj sin 0' + U2 cos 0' 
Uj cos 0' - U2 sin 0’ 


and 


Cl 


CiUi - C4U2 


C3U4 


[vl.vlf^ 


(B 15 ) 


(B 16 ) 


(B 17 ) 


(B 18 ) 
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Similarly, by comparing numerator coefficients. 


(CjCg - C2C4) cos + (c'iCg + C2C3) sin 


, , *,M/2 fhK"l ■ - C 3 U 4 + C2(C3U2 - C’jUj - 

(Uf 1- 

+ [c6(c'iUi - C4U2 - C3U4) + c'lCjC^ sin 4. 


COS ^ 


f C2 COS cp + C5 sin 9 ] cos 9 + C0 sin 9 sin ^ 

^ sin 0' + U2 cos 0') - C2(U3 sin 0' - U4 cos 0') 

^CgUj + C2C2J sin 0' + CgU2 cos 0' 


(uf 


2 u 2'\>/2 


cos 


'2J 


sin 


^C'2 sin 0" - Cg cos 0") cos - Cg cos 0” sin 


02 + U2) 


COS ^ 


+ ^0^2 ^ ^ ^2^3) ^ 


sin ^ 


CjCg - 0304^ sm 4> + (^CjCg + €203^ cos 4> 


H * '^2)'^^ 


Cg(c’iUi - C4U2 - C3U4) + C1C2C3 


cos ^ 


(B19) 


(B20) 


(B21) 
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. jM iff 

Cg sin 0 cos ^ - (^C2 cos 0 + Cg sin 0 ; sin ^ 

= 7 + C2C3) sin 0' + CgUg cos cos # 

(uf + u2)/ 

- |c5(Uj sin 0' + U2 cos 0') - 02(113 sin <p' - U4 cos <p'^ sin (B23) 


-Cg cos 0 cos # - (C2 sm 0 ~ Cg cos 0 ; sm 

= — 777(1^6^2 sin <p' - (CgUi + C2C3) cos <lA cos $ 

(uf . u|)l/2 

- jcg(U2 sin 0' - Uj cos <p') + 02(03 cos 4>' + U4 sin sin ^ (B24) 


Equations (B19) and (B22) can be solved to yield 


r>"r>” r'”n” 

^1^5 " ^2'-4 


u?*u| 


J-^f[c5(ciUi - C^Uj - C^U4) - C2(C\U3 - cIjUj + CiU4)] cos (*’' - i) 


- - C4U2 - C3U4) + CjCjCj] sin (f" - 


(B25) 


^TT_Tf 

CjOg + O2O3 


u? + «i 


Ce(c'jUj - C4U2 - C3U4) + C1C2C3I cos (#" - 4) 


+ 05(o'iU4 - O4U2 - O3U4) - 02(0403 - O3O2 + O4O4) sin (4>" - $) 


(B26) 
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Similarly, equations (B20) and (B23) yield 


f if • ^ff 

P COS 9 + Cp- sin 0 


1/2 I (^5^1 “ ^2^3) (^5^2 + ^2^4) 

(Uf + ui) 


cos ($” - 4>) 


^C0Ui + C2C3) sin <p' + C0U2 cos Ip' 


sin (4“" - 4>) 


and 


(B27) 


_• ^ff 

C0 Sin (p = 


(uf + IL 


[^6^1 ^2^3) ■'■ ^6^2 


cos (<i>" - <i>) 


(C5U1 - C2U3) sin p' + (C5U2 + C2U4) cos p' sin (<i>" - 4>)j 

(B28) 


and equations (B21) and (B24) yield 


^ • JL ^ t ^ t t 

P Sin 9 - Cc cos 9 


rr^T^ F5U2 + C2U4) sin <)■■ - (C5U, - C2U3) cos <!>' 

Pi* 4 ) ^ 

- jCgU2 sin 0' - (CgUj + ^2^3) 


COS (4>" - 4>) 

sin (4»" - 


and 


(B29) 


A.II 1 

-C0 cos P = 


(uf + u|) 


1/2 1 L 


C0U2 sin p' - (C0Ui + C2C3) cos p' cos ("i*" - 4>) 


[(C5U2 + C2U4) sin P' - (C5U1 - C2U3) cos 0*] sin (4>" - 

(B30) 
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When equations (B15) and (B16) and equations (B28) and (B30) are expressed as 
an equivalent ratio, the angle (^" - may be found to be 


($" - $) = tan"l 


C3U2 


U1U4 + U2U3 
Equations (B28) and (B30) also yield 


(B31) 


Cg = * ^2^3) 

(uf.ul)''"’" 


(C5U2 + C2U4) sin (<i>" - 4>) + CgU2 cos (4>” - 4>) 


sin (6" - 6 ’) 


(B32) 


Equations (B27) and (B29) yield 


C2 = 


(uf + u|) 


l'2Vb 


(C5U2 + C2U4) cos (<(>" - 4>) - CgU2 sin (4>" - 4>) 


cos (6” - <t>') 


- ^(CgUj - C2U3) cos (4>" - <i>) - (CgU^ + C2C3) sin - 4>) sin {d>" - d)')| (B33) 


and 


C5 = 


(»5 * «1)‘ " 


(C5U2 + C2U4) COS (‘f>” - 4>) - CgU2 sin sin ((t)'" - 6 ') 

+ jjCgUj - C2U3) cos (<{>" - <i>) - (CgUj^ + C2C3) sin ($” - <i>)| cos (*" - 6' 

and finally equations (B25) and (B26) yield 


(B34) 


Cgfuf + u| 


^4 = . . + C3U4)[c5 cos (#” - $) - Cg sin (#" - <!>)_ 

- C'lRcgUi - C2U3) cos (#" - - (C0Uj + C2C3) sin ($" - 4>) 


+ CjCg^Uj + - C2(C'3U2 - C4U4) cos (4>" - (B35) 
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and 


^=3 = * ^'1/2 {c'l[(C5^Jl - C2U3) Sin (*" - *) * (CeUi + CjCj) con (1." - *)] 




- (C4U2 + C3U4)[cg cos («>" - ^) + Cg sin ($” - f>)] 

- C2(C3U2 - C^U4) sin ($” - 4 >) - C'^C^’ (uf + 


(B36) 


Since equations (B7) and (B8) represent the sequential application of two general 
perspective projection relationships and are found to exhibit the same form of equations 
as a single perspective projection, two (or more) sequentially performed perspective 
projections must be equivalent to a single perspective projection. 
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TABLE L- TARGET POSITIONS DETERMINED FROM DIRECT APPLICATION OF TRANSFORMATION 

TO PHOTOGRAPHS FROM DIFFERENT CAMERAS 


Camera 

format 

size 


Target A 


Target B 


Lateral position, X, m Longitudinal position, Y, m Lateral position, X, m Longitudinal position, Y, mj 


Actual 


Determined 


Actual 


Determined Actual 1 Determined 


Actual 


Determined 


35 mm 

1 

6 cm X 8 cm 
10 cm X 13 cm 
20 cm X 25 cm 


0.920 


0.923 

.915 

.926 

.922 


3.658 


3.652 

3.650 

3.655 

3.662 


1.530 


1.531 

1.532 
1.537 
1.530 


0.914 


0.911 

.913 

.918 

.916 


CO 

to 



TABLE n.- MODEL POSITION^ AND HEADING ANGLE DETERMINED FROM TRANSFORMATION 
ADJUSTED FOR ELEVATION OF MODEL TARGETS 


(a) Far-left model position 


Case 

Lateral position, X, m 

Longitudinal position, Y, m 

Heading angle, i//, deg 

Actual 

Determined 

Actual 

Determined 

Actual 

Determined 

1 

-0.610 

-0.600 

13.411 

13.373 

-30.00 

-30.12 

2 



-.614 



13.375 

-15.00 

-15.60 

3 



-.577 



13.360 

.00 

4.34 

4 



-.601 



13.378 

15.00 

17.62 

5 



-.631 1 



13.446 

30.00 

31.80 


^Model positions are depicted in figure 5, 


(b) Far -central model position 


Case 

Lateral position, X, m 

Longitudinal position, Y, m 

Heading angle, deg 

Actual 

Determined 

Actual 

Determined 

Actual 

Determined 

6 

0.610 

0.600 

13.411 

13.405 

-30.00 

-29.95 

7 



.607 



13.387 

-15.00 

-14.22 

8 



.600 



13.404 

.00 

.60 

9 



.580 



13.406 

15.00 

15.26 

10 



,608 



13.411 

30.00 

30.31 


(c) Far-right model position 


Case 

Lateral position, X, m 

Longitudinal position, Y, m 

Heading angle, xl/, deg 

Actual 

Determined 

Actual 

Determined 

Actual 

Determined 

11 

1.829 

1.765 

13.411 

13.354 

-30.00 

-32.49 

12 



1.829 



13.377 

-15.00 

-12.69 

13 



1.830 



13.398 

.00 

-.42 

14 



1.856 



13.453 

15.00 

15.34 

15 

' 


1.783 



13.344 

30.00 

32.17 
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TABLE II.- Continued 




(d) Intermediate-left model position 


Case 

Lateral position, X, m 

Longitudinal position, Y, m 

Heading angle, i//, deg 

Actual 

Determined 

Actual 

Determined 

Actual 

Determined 

16 

-0.610 

-0.599 

10.973 

10.957 

-30.00 

-30.33 

17 



-.598 



10.968 

-15.00 

-14.26 

18 



-.609 



10.967 

.00 

3.54 

19 



-.616 



10.973 

15.00 

17.67 

20 



-.588 



10.963 

30.00 

32.19 


(e) Inter mediate -central model position 


Case 

Lateral position, X, m 

Longitudinal position, Y, m 

Heading angle, i//, deg 

Actual 

Determined 

Actual 

Determined 

Actual 

Determined 

21 

0.610 

0.588 

10.973 

10.964 

-30.00 

-29.00 

22 



.604 



10.949 

-15.00 

-13.09 

23 



.599 



10.972 

o 

o 

.62 

24 



.604 



10.963 

15.00 

15.99 

25 



.612 



10.974 

30.00 

31.28 


(f) Intermediate-right model position 



Lateral position, X, m 

Longitudinal position, Y, m 

Heading angle, xl/, deg 

Case 








“ 


Actual 

Determined 

Actual 

Determined 

Actual 

Determined 

26 

1.829 

1.806 

10.973 

10.946 

-30.00 

-29.95 

27 



1.800 



10.941 

-15.00 

-14.92 

28 



1.828 



10.951 

.00 

.67 

29 



1.824 1 



10.961 

15.00 

15.62 

30 

1 


1.821 



10.947 

30.00 

29.84 
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TABLE II. - Concluded 


(g) Near -left model position 


Case 

Lateral position, X, m 

Longitudinal position, Y, m 

Heading angle, xj/, deg 

Actual 

Determined 

Actual 

Determined 

Actual 

Determined 

31 

-0.610 

-0.609 

8.534 

8.528 

-30.00 

-31.02 

32 



-.610 



8.534 

-15.00 

-15.76 

33 



-.621 



8.547 

.00 

-.47 

34 



-.591 



8.527 

15.00 

15.68 

35 



-.584 



8.539 

30.00 

29.81 


(h) Near-central model position 


Case 

Lateral position, X, m 

Longitudinal position, Y, m 

Heading angle, xp, deg 

Actual 

Determined 

Actual 

1 

Determined 

Actual 

Determined 

36 

0.610 

0.603 

8.534 

8.530 

-30.00 

-30.26 

37 



.612 



8.524 

-15.00 

-16.66 

38 



.606 



8.533 

.00 

.38 

39 



.610 



8.533 

15.00 

15.19 

40 



.615 



8.518 

30.00 

29.46 


(i) Near-right model position 


Case 

Lateral position, X, m 

Longitudinal position, Y, m 

Heading angle, i//, deg 

Actual 

Determined 

Actual 

Determined 

Actual 

Determined 

41 

1.829 

1.814 

8.534 

8.535 

-30.00 

-29.11 

42 



1.818 



8.531 

-15.00 

-14.00 

43 



1.830 



8.533 

.00 

.86 

44 



1.811 



8.520 

15.00 

15.85 

45 



1.817 



8.522 

30.00 

30.88 
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Figure 1.- Tilted perspective 


optical axis 


Principal plane 



Lane 


n with coordinate rotations. 






(b) Image plane. 

Figure 3.- Views of reference grid in object and image planes. 
Rotated principal axes are denoted by bars. 


4 ^ 

Oi 


Figure 4,- View of grid and targets for direct application of transformation. Dimensions are given in meters. 








NASA- Langley, 1976 L-10411 






(d) Intermediate -left position. 




(f) Intermediate-right position. 



(g) Near -left position. 



(h) Near -central position. 



(i) Near -right position. 


L-76-121 

Figure 5.- Positions of landing gear model with raised targets. Model is shown with a 30° clockwise 
heading {\{/ = -30®). Corners of reference grid are indicated by circles and the specified grid 
coordinates are given in meters. 
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